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Abstract 

The problem of phase space transport which is of interest both theoretically and from the point of 
view of applications has been investigated extensively using geometric and probabilistic methods. Two 
of the important tools for this that emerged in the last decade are the Finite time Lyapunov exponents 
(FTLE) and the Perron-Frobenius operator. The relationship between these approaches has not been 
clearly understood so far. In this paper a methodology is presented to compute the FTLE from the 
Perron-Frobenius operator, thus providing a step towards combining both the methods into a common 
framework. 

Transport on finite time scales in dynamical systems has been studied intensively in the last 
two decades. A significant motiviation for this study has been transport in fluid systems vary- 
ing from micro fluidics to geophysical flows. These studies received a significant boost with the 
development of the concepts of finite time Lyapunov exponents (FTLE) and Lagrangian coher- 
ent structures (LCS), p], [2], [3], which act as transport barriers. A probablistic approach to 
transport in dynamical systems developed in the last decade, which uses the Perron-Frobenius 
(PF) operator, pQ, to identify the so called almost invariant sets, [5j, |6j. The relationship be- 
tween the two methods has not been well understood though both methods often yield roughly 
the same results, [7]. In this paper a method to compute the FTLE from the PF operator 
is given by using an alternative definition of the FTLE based on the covariance of probabil- 
ity density functions. This method of computing the FTLE from the PF operator eliminates 
the need for long time integrations of trajectories, often a source of errors. The method is 
illustrated with two simple examples of two dimensional fluid systems. The redefinition of 
the FTLE along with its computation from the PF operator are a step towards combining the 
geometric and probabilistic methods into a common framework. 
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1 Introduction 



The problem of phase space transport has important applications such as in mixing and separation problems 
in fluid flows that vary in scale from the micro to the geophysical, interplanetary transport and instability of 
mechanical systems, to name a few. A variety of dynamical systems methods have been studied over the past 
three decades to explain transport mechanisms, to detect barriers to transport, and to quantify transport 
rates, see eg pQ, 0, 0, [S], [S], [§], [ID], [H], [H], [IS]- These methods fall into two main categories, the 
geometric and the probabilistic. Under the umbrella of geometric methods are the techniques of invariant 
manifolds (of fixed points) , lobe dynamics and Finite time Lyapunov exponents (FTLE) and Lagrangian 
coherent structures (LCS). The method of FTLE and LCS has proven to be particularly useful in studying 
transport in time dependent systems and has found a variety of applications, for e.g., [14] . [15] . |16j . |17j . 
[15] and [IH] ■ The probabilistic approach studies the transport of densities and the so called almost invariant 
sets and coherent sets. These methods too have been successfully applied in the study of various geophysical 
flow problems, |12) and mixing in micro channels [20] . 

The method of LCS studies stretching and contraction around reference trajectories. The LCS method 
is therefore local in nature; it provides information about invariant manifolds that determine transport in 
phase space. The current method of LCS relies on computing the FTLE field using long time computations, 
since LCS usually can be identified only after a significant time of integration. The disadvantage of this 
is that, excessive stretching of line elements can introduce computational errors. Moreover specific checks 
on whether the stretching of line elements is within the linearized regime can be difficult to incorporate in 
the current algorithms on computing the FTLE. The probabilistic method on the other hand ignores the 
local transport structures, but using the transfer operator divides the phase space into maximally invariant 
sets. There have been a few attempts to explore the relationship between the geometric and probabilistic 
descriptions of phase space transport, such as [7] and |21| . The aim of this short paper is to present a 
technique of computing the FTLE using the Perron-Frobenius operator that is a step towards combining 
the geometric and probabilistic methods, by making the Perron-Frobenius operator the common tool to 
both. By utilizing the Perron-Frobenius operator to compute the FTLE, this method also strengthens the 
probabilistic interpretation of the FTLE, identified in [21] and (22[. This approach has the added benefits of 
eliminating the issue of linearization around a reference trajectory in the existing formulation of the FTLE 
and reduces the time to compute the FTLE field for time-dependent and periodic flows. 
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2 Review of FTLE and LCS 

The formulation of FTLE and the Pcrron-Frobenius operator is reviewed in this section. This review is 
intended to provide an intuitive background and set the context for the computation of the FTLE using the 
Perron-Frobcnius operator. For the rigorous definitions and details on these methods, the reader is referred 

to0, ©, 0, 0, ffl. 

2.1 Finite Time Lyapunov exponents 

Let M C W 1 be compact and §\ a (x) :MxRxRi->Mbea smooth flow on M. Let the associated vector 
field be x = F(x,t) for x g M. Consider a reference trajectory passing through the point x(to) and a 
perturbed trajectory passing through the point x(to) + Sx(to) at time to- The flow (b\ Q maps these points to 
4>l Q (x(to)) and 0' o (x(£o) + Sx(to)) at time t and the perturbation grows to Sx(t). 
Expanding </>j o (x(£o) + Sx(to)) in a Taylor series about the point x(£q) we get 



6x(t) = <4( X (M) - <(x(t ) + 6x(t )) = -^«x(to) + O(||<5x(t )|| 2 ) 



(1) 



The norm or magnitude of 5x(to + t) can be found using the standard inner product on 
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(2) 



where * denotes the transpose and the gradient I 1 is evaluated at x(to). The maximum growth of a 
perturbation is therefore given by the maximum principal stretch, i.e., by the maximum eigenvalue of C. 



max \\Sx(t)\\ = y/X max (C(x(t ),t ,t) \\Sx(t )\\ 



(3) 



where C(x(to), to, t) ~ (^^j (^ET J is the Cauchy-Green tensor. The growthin the peri urbation depends 
on the initial point x, initial time to and the evolution or integration time T — t — to- 



Definition 2.1 The maximum FTLE is defined as, Jlf, f3jj, 



a(x(t ), t , T)=^ log ^X max (C{x(to),to,t)) 



(4) 



The leading FTLE gives the time averaged rate of linearized stretching in a neighborhood around a reference 
trajectory. It is intuitively clear that regions of the phase space with locally high values of FTLE will stretch 
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and separate. The sets with high FTLE act as repelling barriers in the flow. This intuitive idea of barriers is 
formalized by the concept of Lagrangian coherent structures (LCS), defined as ridges in scalar FTLE field, 
[2], [3j. Ridges can be defined precisely by appealing to differential geometric quantities as in [23], [2] and 

2.2 Perron-Frobenius operator 

Let \x denote the Lebesgue measure and let B be a measurable set and / <G L 1 be a probability density 
function, L 1 being the space of Lebesgue measurable functions. The unique operator V\ : L 1 n- L 1 defined 

by 

/ V\jd i i= [ fdn (5) 

is called the Perron-Frobenius operator for the flow (j)\ o , [4j . Equation [H] which holds for all /i measurable 
sets follows from the Radon-Nykodym theorem. 

In practice it is usually necessary to numerically approximate the operator V\ . This is done by discretiz- 
ing the domain, M, into a finite number of sets, say {B\, £?2, B n } which is essentially a grid of boxes. A 
projection ir : L 1 i— > span{Bi, ...,B n } defined by nf = X)"=i C i^-Ba where Xb ( is the characteristic function 
of the set Bi and Ci — gives a finite dimensional approximation of /. Since / is a probability density 

function Cj = ^ B .y Similarly V\ a f is projected on span{Bi, ...,B n }. The operator P t * o : nf H> ^V\ a f is a 
linear operator between finite dimensional vector spaces. Further taking the box measures [i{Bi) = p{Bj) 
for all i,j € [l,n], P* becomes a stochastic transition matrix. The entries of the matrix P are determined 
by a Monte-Carlo simulation [5j and . Each box in the domain contains a fixed number of points (initial 
conditions), which are integrated from a time to to t. The final position of the points gives the matrix P as - 



A time reversible operator P is required to apply the above definition for flows in forward time, [6]. This is 
achieved by creating a reversible Markov operator P r given by 



where P is the time reversed analogue of P. Its elements are given by 



(W )ij - U} ^ hl (8) 
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where uij and Uu are components of the first left eigenvector, m, of P[ ' . For a volume preserving flow in 
which the domain is uniformly discretized, Pf = (P/ )*, the transpose of P* . Henceforth for convenience the 
time reversible operator (P/ ) r is referred to simply as P unless there is an ambiguity about the initial and 
final times to and t. The Markov operator P has the semigroup property of P\ Q — P t s Pj, where s € (to,t). 
This property is used to simplify the computations of the FTLE. 



3 Computation of FTLE using the Perron- Frobenius operator 
3.1 Definition of set oriented FTLE 

The concepts of FTLE and LCS reviewed previously have been used fruitfully in many areas as has been 
pointed out earlier. However the standard computational implementation of FTLE, using finite differences, 
suffers from some drawbacks, stemming from the integration time T and the linearization around reference 
trajectories. In equation [l] it is assumed that the second (and higher ) order terms are negligible. However 
the magnitude of the higher order terms depends on the evolution time T = t — tg. If the evolution time 
T is too high then the higher order terms may be comparable to the first order terms in equation [l] If the 
evolution time T is too low, then one cannot detect any interesting structure in the FTLE field. Often, 
the evolution time T is selected in a subjective fashion without checking the validity of equation [T] The 
problem with this is illustrated with the example of the double gyre flow, a prototype in the LCS literature, 
[3J [7]. The double gyre flow is defined by the stream function ip(x,y,t) — A sin(7r/(x, t)) sin (iry), with 
f(x, t) = ex 2 smtut + x(l — 2e sinwi), where we use the parameters A = 0.25, u> = 2ir and e = 0.25. The time 
period of the flow is t — 1. Now consider the equation [T] with the second order terms explicitly written, 

«x(*o +t) = <(x) - <(x + <5x) = ^x(io) + ^x(t )*^«5x(t ) + O(||<5x(i )|| 3 ). (9) 
In this equation it can be shown that the second and first order terms are comparable in magnitude. The 



shows the plot of A max I I ||<5x|| and 



max-norm for matrices is used for this comparison. Figure 1 

|V„ (^) |fa|>. The m .^t„de of the second order ter ms is more than half that of the first order 
terms in equation [9] This is true for a smaller evolution time as well. The ridges in FTLE field for the 
double gyre flow are generated by moving instantaneous stagnations points (ISPs) . Material line elements 
close to the ISPs stretch in a shorter time as compared to material line elements farther away. Therefore 
ineresting structures such as ridges in the FTLE field require a long period of integration during which the 
nonlinear deformation of line elements close to the ISPs occurs. 
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Figure 1: (a) Comparison of first- and second-order terms in eq. q9p for T = t — to = 10. 

To study the linearized flow in equation [9j numerics have to be carefully planned to ensure that the 
approximation of linearization is valid. Usually this is accomplished by selecting a very small 5x and either 
by keeping the time of integration T small enough or rescaling the perturbation as it grows very large. 
If perturbations around a specific trajectory grow too large and need to rescaled, then additional reference 
points have to be introduced in that region to obtain the FTLE field at a fine enough resolution. Alternatively 
the computations could start with a crude mesh of initial points and refined iteratively by introducing new 
initial points based on the finite time stretching, requiring adaptive meshing of initial conditions, which was 
explored in [17j . However such techniques of mesh refinement are not based on the magnitude of the second 
order terms in equation [9j Beyond the fact that numerical algorithms for the calculation of FTLE do not 
check the validity of the linearization of the flow, there is the additional fact that interesting structures in 
the FTLE field develop for longer eveolution times during which material lines stretch nonlinearly and even 
fold as shown in figure [2] which was discussed earlier in |17j . 

The difficulties in the computation of the FTLE can be overcome with a set oriented definition of the 
FTLE. Instead of tracking only a reference trajectory and the adjacent nodes of the finite difference mesh, 
one tracks the movement of the whole set (shown in gray) in figure [2] and consider the distribution of points 
in the set, then folding and nonlinear deformation can be accounted for. A new definition of the FTLE 
was proposed in [3T] which computed the deformation of sets instead of the stretching of line elements and 
partially addressed the concerns of second order terms, integration time and folding of material lines. We 
review the new definition here. 

For illustrating the concept we assume the flow <j>\ : M M is is two dimensional with M C M 2 . The 
method of computing FTLE using the SVD of the Cauchy-Green tensor essentially computes the linearized 
stretching of a neighborhood, a set B, under the action of the flow as shown in figure [3j The FTLE for 
the reference trajectory in this case is a = i log (—J where T is the time of evolution of the trajectory. 

In the standard FTLE method the evolution of the set B is assumed to be determined by the evolution 



G 




Figure 2: Folding of material elements. 




B 

Figure 3 : Deformation of a blob under the flow 

of two vectors, the directions of principal stretches. For the set oriented definition the deformation of the 
set B is tracked by the evolution of a random vector X = [XijX^] defined by a probability density function 
f(xi,X2) which is initially a uniform probability density function supported on B given by, / = }„. Xg, 
where Xg is characteristic function of B. The covariance matrix of / is Zy = E[(Xi — Xi)(Xj — Xj)}, 
with i = l,2 and j = 1, 2 where [Xi, X 2 ] is the mean value of the random vector X and E[-] denotes the 
expectation with E[X] = [J X\f{xx, X2)d/j,, J X2f(xi,X2)dfi]. Under the action of the flow (jr tQ , f is mapped 
to V\ f where V\ is the associated Perron- Frobenius operator. In matrix notation the covariance Iij is 



E[(X 1 -X 1 )(X 1 -X 1 )] E[{Xi-X 1 ){X 2 -X 2 \ (|(() 
E[(X 2 - X 2 )(X 1 - X{)] E[(X 2 - X 2 )(X 2 - X 2 ) 

Definition 3.1 Let 1(f) be the covariance of f andl(Pf) the covariance of Pf and let \ max (I) denote the 
maximum eigenvalue of I . Then the FTLE of B denoted by ai(B,to,t) is defined as - 



ai(B,to,t) = — log f V ^gggSB S ) (ii) 
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It can be shown that by direct calculation that the covariance FTLE obtained from Definition 3.1 and the 
standard FTLE have the same value if the second and higher order terms in equation ^ are negligible when 
compared to the first order terms. In this case an initial circular blob deforms into an ellipse as shown in 
figure [3J The eigenvalues of the C are a\ and a\ . The standard FTLE is a — -j^-j- log ^ . The covariance 
matrix for the deformed ellipse is 



I{Pf) 



( 1 , \ 

j7raja 2 

\ \-na x a\ J 



(12) 



giving X ma ,x(I(Pf)) = 1 [ 7ra i a 2- The covariance matrix /(/) is a diagonal matrix with 7(/)n = — j?™ 4 
giving X m ax(I{f)) — jTra 4 - This gives the covariance FTLE, 07 = j^-j- log f v^i° 2 ^ p or volume preserving 
flows, -Ka\a,2 — Tra 2 . This gives the covariance FTLE, a = j^- log the same value as the standard FTLE. 
The standard FTLE is obtained by the stretching of the value of the covariance based FTLE is the same as 
that obtained from traditional FTLE calculation using line stretching. 

Since I is the covariance of /, it provides a probabilistic interpretation of the FTLE. Simultaneously 
I can also be interpreted as the moment of inertia of the set B and provides a geometric description of 
the deformation or distortion of the set. The definition of 07 avoids the linearization of the flow and the 
computation of the stretching of line elements. Further it is a set-oriented method and directly computes 
the deformation of a set instead of inferring it from the deformation of line elements. If the linearized 
flow is valid the covariance FTLE and the standard FTLE are equal. However the covariance FTLE is a 
better measure of finite deformation of sets. The covariance based method of FTLE is a bridge between the 
geometric approach of measuring line stretching and the probabilistic approach of almost invariant sets. 

In practice any set B and a density function / supported on B has to be approximated by discrete 
points. If the function / = jjgs %b is & uniform density function, then B can be approximated by N points, 
with coordinates (xu,X2i), i ranging from 1 to TV. Each of these N points has discrete measure jy^g) ■ 
Similarly <j)\ (B) is also approximated by the images of the iV points, (cfil (xu,X2i)), each with measure 
jv>(<4 (B)) • ^ 0r vomme preserving flows this is the same as jy^uy- The covariance for the function Pf, 
I(Pf), is approximated by 



I(Pf) 



Eti($o(*ii) -£i)) 2 TtMMi) -zi)($ (z 2i ) -l 2 ) 



N (At 



(13) 



where £ 2 ] is the mean or expected value of the discrete approximation of (f>* (B). The initial covariance 
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matrix can also be obtained similarly. The maximum eigenvalues of /(/) and I(Pf) can be plugged into 
equation (3.1 1 to give the covariance FTLE. 



3.2 Computation of covariance FTLE using Perron-Frobenius operator 

The covariance FTLE is a set oriented redefinition of the FTLE. Therefore is natural to expect the Perron- 
Frobenius operator play a role in the computation of the covariance FTLE. This will be particularly useful 
since only short time integrations are necessary to calculate the Perron-Frobenius operator P. The operator 
P t * can be found using a suitable set of short time intervals, [toj^i]; [^i , ^2] , [inj*]; Pt = PuPtl—Pt n - Let 
/ be a uniform density function supported on a set B. Then as before Pf maps / under the action of the 
flow to Pt f, P' '■ f l— ► Pt f- For time independent flows or periodic flows with period At, this becomes 
particularly easy, with P t * o = (p*°+ A *)«+ 1 where At = t i+ i — ti for i = 0, l,...,n. This avoids long time 
integration of trajectories. 

For two dimensional flows, (j>\ Q : M 2 H> M 2 , the operator V\ is numerically approximated by P* using 



the box discretization method described in section 2.2 Each box contains a fixed number of uniformly 
distributed points. The entries of the matrix P t ' o 0+At are calculated or a suitably chosen At and the time 
reversible operator is found from equations ([6]), Q and ([8]). To compute the FTLE over the entire domain, 
we take the projection of a uniform density function such that nf = u % = [0 ... l//i(Pi) 0], with the ith 
column being equal to the inverse of the measure of Bi. The vector u l evolves to v l (to + At) at to + At, 
given by, 

v*(t + At) = u* P*°+ A * = -L-(P^% (14) 
Similarly the evolution of Ui at time t = to + (n. + I) At is given by 

v\t) = u* (P^ +At ) n+1 = (P!: +A X +1 (15) 



Equation [15] defines how a set of points uniformly distributed over the set Bi at to is disperese at time t. 

t 



The the covariance of this set can be found from the ith row of (P t ' 0+A *) n+1 as follows 



I zf^i-tjiti-f^ptx) EUaci-tnpioh)- 

where [Ci, C2] are tne centers (mean values) of each of the boxes, Bj and [£1,^2] = SjLi([Ci) C2\(Pt° +At )ij)- 
Intuitively the uniformly distributed points in box Bi are mapped into some of the N boxes by the flow 
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Bi 



► 




Figure 4: Illustration of the box discretization method. BoxB; with measure is mapped into three boxes Bj , and B m by 
the flow map </>(£?)?. The values of (P} )ij, (Pt )ik an d (Pt )im are given by equations lol [7l and fsl 




Figure 5: Illustration of the calculation of the covariance, I(v l (t)) in equation |l6| 

map 4>l , say boxes Bj , Bj~ and B m as shown in figure [4| If the number of points that are mapped into each of 
the three boxes are Nj , Nk and N m , then each of the three sets Bj-n0| o (Bj), BkD(f>l (Bi) and B m n<f>l (Bi) are 
approximated by Nj , Nk and N m uniformly distributed points in the boxes Bj , Bk and B m . The mean or av- 
erage of these dispersed points is ^u^]. The vector v l (t) = ^y[0, (Pt o )ij,0, (P&ik, 0, (-P'Jim, ■■■0] 
with the nonzero values in columns j, k and m. The mean or center of this set is given by ^ J 2 ]- Proceeding 
thus, the covariance FTLE for each of the boxes in the domain can be found from equation |l6]). 

It should be emphasized that the computation of the FTLE using the Perron-Frobenius operator requires 
a set oriented definition of the FTLE. One cannot use the operator P to track the evolution of sets of zero 
measure, since the integral on the right hand side of equation [5] is zero for this case. The standard method 
of FTLE which requires the evolution of individual trajectories passing through the nodal points of a mesh 
at the initial instant of time, does not make any direct reference to sets of positive measure around these 



10 



nodal points. Hence a set oriented definition of FTLE is necessary to take advantage of the Perron- Frobenius 
operator. In practice the computational approach for this new definition too has to make use of discrete 
initial conditions and trajectories , but sets of positive measure are explicitly modeled by these discrete 
initial conditions. 

4 Examples 

In this section we illustrate the computation of the FTLE using the Perron-Frobenius operator for two flows. 
The first is the lid-driven cavity flow studied for its mixing properties in [20] and the double gyre flow that 
has been a prototype flow in the LCS and almost invariant sets literature, [5] , [T7] , [7] . 

4.1 Lid driven cavity flow 

The problem of transport in the model of the lid driven cavity flow has been investigated in |20j and 
is considered here as the first example because of its simple piecewise steady velocity field. The flow is 
described by the stream function 

2 

V^,y) = 5>„C„/„(y)sin(^) (17) 

n— 1 

defined on the domain [0, a] x [—6, 6] for time > t > Tt/2. where 



and 




For time T//2 > t > Tf, the sign of the velocity term XJ\ is changed. This reflects the streamlines about 
x = a after a time Tf/2. 

Using symmetry arguments given in [23] and |25) . a specific ratio of the magnitudes of the terms U2/U1 
along with a fixed value of the period of the flow Tf is found such that it generates three period-3 fixed points 
in the domain [0,a] x [—6,6]. The specific values of the constants that we borrowed from [24] and [25] are 
Ui = -1 and U 2 = 0.841298 and t} /2 sa 4.740202 for the domain [0, a] x [-6, 6] = [0, 6] x [-1, 1]. 

A perturbation of the time period of the flow from the critical value of tJ/2 sa 4.740202 destroys the 
fixed points. The method of almost invariant sets was used in [26j and [20j to study mixing for different 
values of the perturbed time period. Since our main interest here is to illustrate the FTLE computed from 
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the Perron- Frobenius operator, we choose a single case, a specific value of half time period r//2 = 4.848 for 
our study. 

The standard FTLE computation was done for the system for different integration times; Tf, 2rf and 
3r/, the results of which are shown in figure |6ja)-(c) by integrating points initially spaced at a distance of 
0.05. The domain is then divided into 4800 boxes each box containing a 100 points and the time reversible 
Perron-Frobenius matrix P r is found by integrating a total of 480,000 initial points for a period equal to 
Tf. The vector space has 120 x 40 dimensions which are the discretized finite approximation of the infinite 
dimensional space of Lebesgue integrable functions. Uniform density functions each supported on a box 
are the bases vectors for the finite dimensional vector space. The covariance based FTLE computation was 
performed using the Perron-Frobenius operator Pq f . Figure [6] (d) is obtained by computing the covariance 
of the bases functions from the rows Pq 1 , while figure [6] (e) and (f ) are obtained by computing the covariance 
of the basis functions from the rows of Pq T/ = {Pq S ) 2 and Pq T/ = (Pq 1 ) 3 respectively. 

It is evident from the figure [6] that the FTLE field has the same ridge structure when computed by the 
line stretching approach or the covariance approach using the Perron-Frobenius operator. Computing the 
covariance based FTLE field directly by integrating all the initial conditions to 3r, while more accurate is 
computationally more intensive. In fact the time taken for this method increases almost linearly with the 
integration time. The FTLE field computed from the Perron-Frobenius operator has the same features as 
the covariance FTLE field by integrating the points in each of the boxes for 3r and computing the covariance 
of the bases functions which is shown in figure [7] 

While the FTLE field computed directly by integration and that computed by the Perron-Frobenius 
operator have the same ridge features, they differ in the magnitude of the FTLE field. This discrepancy in 
the magnitude is due to two reasons. The Ulam method of approximating the Perron-Frobenius operator by 
a matrix introduces approximations. The measure of each of the boxes is approximated by a finite number 
of discrete points introducing further errors. Lastly the calculation of the covariance matrix I(v l (t)) using 
equation as explained in figure [4] introduces another level of approximation. As the box size in the compu- 
tation is reduced and/or the number of points per box is increased, it is natural to expect that the matrix 
approximation of P and the covariance computation become more accurate, though no rigorous proof is 
offered here. 

4.2 Double gyre flow 

The double gyre flow which is time dependent but periodic, has been one of the prototype flows in the LCS 
literature, [2], [7]. The double gyre flow is defined by the stream function ip(x, y,t) — A sin(7r/(x, t)) sin (ny), 
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in 



(a) T = Tf 



(d) T = T, 





(e) T = 2r f 




(c) T = 3r f 



(f) T = 3r f 



Fi gllXG 6i FTLE for integration time T for the lid driven cavity flow. On the left (a-c) is the FTLE field calculated using the 
standard method of stretching of line elements. On the right (d-f ) is the FTLE calculated using the operator P and equation 

Jlet. 



with the parameters A — 0.25, ui = 0.2ir and e = 0.25. The time period of the flow is r = 1. The domain 
is discretized into 50000 square boxes of size 0.005 each containing 625 points. The FTLE field computed 
using the stretching of line elements and the FTLE field computed from the Perron-Frobenius operator is 
shown in figure [8] 

As observed in the discussion on the lid driven cavity flow, the FTLE field has the same ridge features, 
when the covariance of the sets Pf are computed directly by integrating all the points in each box to T = 10, 
as shown in figure [9] and when computed using the Perron-Frobenius operator. Moreover because of the more 
accurate approximation of the domain by a higher number of boxes, the covariance FTLE field computed 
by a direct integration and that computed by the Perron-Frobenius operator are closer to each other in 
magnitude. 
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2 4 6 



(c) T = 3t/ 

Figure 7: Covariance FTLE for lid driven cavity flow computed by direct integration of initial conditions using equation | |13| 
for various integration times. 




1 2 



(c) T = 10 (d) T = 10 

Figure 8: FTLE for integration time T for the double gyre flow. On the left; (a) and (c) is the FTLE calculated using the 
standard method of stretching of line elements. On the right; (b) and (d) is the FTLE calculated using the operator P and 
equation ( |16| 
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Figure 9: Covariance FTLE obtained by direct integration for an integration time T = 10 using equation ( |13| for the double 
gyre flow. 

5 Conclusion and Discussion 

A method to compute the FTLE field using the Perron Frobenius operator has been introduced in this 
paper. To do this a modified definition of the FTLE was used which identifies deformation of sets of positive 
measure instead of the stretching of line elements. This methodology is general enough to be used in time- 
dependent flows and is particularly useful in the case of time-independent and periodic flows, where it offers 
a significant computational advantage by eliminating long time integrations. The method of computing the 



covariance FTLE from Perron- Frobenius operator using equation ( 16 ) was illustrated by two examples; the 
lid driven cavity flow and the double gyre flow. In both cases it was shown that the FTLE field computed by 
the new method has the same topological features as the FTLE field computed using the standard approach 
of line stretching and the direct computation of the covariance FTLE from the dispersion of points as in 



equation (13 1. The computational time required for the calculation of the covariance FTLE directly from 
the dispersion of points scales linearly with time in the most ideal case. So the time taken to calculate 
(Ji{to, to +nAt), for an integration time nAt, is n times that of the time taken for calculating a i (to, to + At). 
Using the PF operator approach the time necessary to calculate the FTLE is almost independent of the 
multiple n. Thus the method proposed in this paper is n times computationally efficient. 

Comparing the method proposed in the paper, with that of the standard FTLE, the computation of the 
standard FTLE is faster. This is because the standard FTLE computations were done a very coarse grid in 
both the examples presented. Even though the time of integration is shortes, fewer initial conditions have 
to be integrated and thus the method is faster. However as has been pointed in the paper in section 3.1 this 
is achieved at the cost of erroneously ignoring the higher order terms in equation ([!]) . The covariance based 
method of computing the FTLE does not use the linearized equations of a flow making it more general in 
scope of application. 

Moreover the covariance FTLE gives the deformation of sets a probabilistic interpretation. The new 
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method of computing the covariance FTLE introduced in this paper is based on the Perron-Frobenius 
operator; a common tool in probabilistic methods of phase space transport. This puts the Perron-Frobenius 
operator at the center of study of phase space transport. This is perhaps intuitively obvious since the Perron- 
Frobenius operator contains all the information of global transport of sets. The method of computing the 
covariance FTLE using the Perron-Frobenius operator is a step towards combining the probabilistic and 
geometric methods of phase space transport, into a common unified framework. 
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